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2 ! Abstract 

We study the effect of a magnetic field on the pairing dynamics in two-flavor color supercon- 

ducting dense quark matter. The study is performed in the weakly coupled regime of QCD at 

asymptotically high density, using the framework of the Schwinger-Dyson equation in the improved 

rainbow approximation. We show that the superconducting gap function develops a directional 

dependence in momentum space. Quasiparticles with momenta perpendicular to the direction of 

, the magnetic field have the largest gaps, while quasiparticles with momenta parallel to the field 

00" 

have the smallest gaps. We argue that the directional dependence is a consequence of a long-range 



interaction in QCD. The quantitative measure of the ellipticity of the gap function is determined 
by a dimensionless ratio, proportional to the square of the magnetic field and inversely proportional 
to the fourth power of the quark chemical potential. For magnetic fields in stars, B < 10 18 G, 
the corresponding ratio is estimated to be less than about 10 -2 , justifying the use of the weak 
magnetic field limit in all stellar applications. 
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I. INTRODUCTION 



Quantum chromodynamics (QCD), the fundamental theory of strong interactions, pre- 
dicts that quark matter at sufficiently high densities and sufficiently low temperatures is 
a color superconductor flfl. (For reviews, see for example Refs. p, |7J.) The property of 
asymptotic freedom in QCD ensures that such matter is weakly interacting at asymptoti- 
cally large densities and, therefore, allows a rigorous treatment of the corresponding non- 
perturbative dynamics of Cooper pairing js-15]. The simplest color superconducting phases 
correspond to spin-zero pairing. Depending on the number of quark flavors participating 
in pairing, one can have a two-flavor color superconducting (2SC) phase % |5| or a color- 

I L 

flavor-locked (CFL) phase [16J. Many additional complications arise when /3-equilibrium 



and neutrality of quark matter is enforced 



17^20). 



Recently, the st udy of color superconductivity in the presence of magnetic fields attracted 



a lot of attention 



21 



lN25|. This interest is primarily driven by potential astrophysical ap- 
plications, where magnetic fields play an important role. In the case of neutron stars, for 
example, the surface magnetic fields can reach up to about B ~ 10 12 G 26]. For magnetars, 
the corresponding fields can be still a few orders of magnitude larger, i.e., B ~ 10 14 — 10 15 G, 
and perhaps even as high as 10 16 G 27J. Furthermore, it is possible that the magnetic fields 



26 



28|. 



in the stellar interiors are much higher and reach up to about B ~ 10 18 G 

In order to understand the properties of two- and three-flavor color superconducting 
phases with spin- zero pairing in a magnetic field, it is important to first recall their electro- 
magnetic properties. Despite being color superconductors, these phases can be penetrated 



by long-range "rotated" magnetic fields, which are not subject to the Meissner effect [291130]. 
The rotated gauge fields are linear combinations of the vacuum photon and one of the gluons. 
While all Cooper pairs are neutral with respect to the corresponding rotated electromag- 
netism, the individual quark quasiparticles carry well defined charges. It is not surprising, 
therefore, that the diquark pairing dynamics is affected by the presence of a magnetic field. 
The recent studies revealed many interesting qualitative features of the magnetic 2SC and 



CFL phases 
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23). However, all such studies share a common shortcoming: they are per- 



formed in the framework of Nambu-Jona-Lasinio (NJL) models with contact interactions. 

In this paper, we extend the analysis of two-flavor color superconductivity in a magnetic 
field by taking into account the long-range interaction in quark matter. In particular, we 
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perform the study in the framework of the Schwinger- Dyson equation for the gap function in 
the weakly coupled regime of QCD at large densities. The long-range interaction is provided 
by the one-gluon exchange, in which the dominant screening and Landau damping effects 
are included. Our study reveals a qualitatively new feature of the magnetic 2SC phase, a 
directional dependence of the gap function, which is a consequence of the nonlocal interaction 
in quark matter. 

In the weak magnetic field limit, we find that the effect of a nonzero field can be mim- 
icked by an effective increase of the strong coupling constant that governs the Cooper pair- 
ing dynamics: g 2 — >■ g 2 (l + esm 2 9 Bk ), where Obu is the angle between the quasiparticle 
momentum and the direction of the magnetic field, and the dimensionless quantity e is a 
measure of ellipticity of the gap function. The latter is given by the dimensionless ratio 
e = 27ir(eB) 2 /(2g 2 p, A ), where B is the magnetic field and p, is the quark chemical potential. 
As one can easily check, this ratio is much less than 1 even for the strongest possible fields 
in stars and, therefore, the use of the weak magnetic field limit is justified for all stellar 
applications. For completeness, we extend our analysis to the case of superstrong magnetic 
fields and find that the value of the gap increases with the field also in this regime. As 
expected on general grounds, the effects of nonlocality of the interaction become negligible 
in superstrong fields and the directional dependence of the gap disappears. It should be 
remarked, however, that our analysis in the case of strong fields is performed with less rigor 
because the gluon screening effects in this case are not well known. 

The rest of the paper is organized as follows. In Sec. [TT], we introduce the model and 
review how a constant rotated magnetic field enters the Lagrangian density and how the 
rotated electric charges of quasiparticles are defined. Explicit expressions for quasiparticle 
propagators in sectors with different rotated charges are presented in Sec. IIHI Then, in 
Sec. IIVt we derive the gap equation and solve it approximately in the limit of a weak 
magnetic field. In the same section, we also obtain an estimate for the gap in the strong 
field limit. In Sec. |Vl we discuss the results and give a brief outlook. Several Appendices at 
the end of the paper contain many technical details and derivations used in the main text. 
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II. MODEL 



As stated in the introduction, the analysis in this study is done in the framework of weakly 
interacting two-flavor QCD at large densities. The quadratic part of the corresponding 
Lagrangian density of quarks in an external rotated magnetic field is given by 



pern 

quarks 



(2.1) 



where is the rotated massless U(l) em gauge field. This field is a linear combination of 
the vacuum photon A, and the eighth gluon G 8 : = cos 6 A ^ — sinOG 8 , where cos 6* = 



g/ \/g 2 + e 2 /3 |29|, |30|. (Here we use the standard convention for SU(3) C color generators 



in the adjoint representation 30].) The quarks carry flavor and color indices ipi ai where 
i G (u,d) = (1,2) is the flavor index and a e (r,g,b) = (1,2,3) is the color index. The 
multicomponent quark spinor field ip is assumed to have the following explicit form: 



( tpur \ 
^ug 
^ub 
^dr 

\ ^db ) 



(2.2) 



Here we assume that up and down quarks have the same masses (m u = rrid = m). In the 
2SC phase, the matrix of chemical potentials £i can have a nontrivial color-flavor structure. 
When /3-equilibrium and neutrality of quark matter is imposed 18), the matrix elements of 
pL read 

2 

/%,a& = [fJ'Sij — fA e (Qf)ij] ^ab + -^^ij{T^)abi (2-3) 

where only one out of three parameters (/i, /x e , and /i 8 ) is truly independent, while the other 
two must be adjusted to achieve color and electric neutrality. For subtleties regarding the 
color neutrality; see Ref. [31 1 . 

The explicit form of the quasiparticle charge operator Q, that corresponds to U(l) era 
gauge group, is given by Q = Qf <g) I c — If <g> (^)c, where Qf = diag(|, — |) is the usual 
matrix of electromagnetic charges of quarks in flavor space, and T 8 is the eighth generator 



of the SU (3) c gauge group in the adjoint representation. In units of e = eg/ g 2 + e 2 /3, the 
Q charges of quarks are given in Table [H 
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TABLE I: Q charges of quarks measured in units of e = eg/y/ g 2 + e 2 /3. 



u r 


Ug 


u h 


d r 


dg 




+2- 


+ h 


1 


i 

2 


1 

2 






In order to simplify the explicit form of the quark propagators in the magnetic 2SC 
phase, it is convenient to introduce the following set of projectors onto the subspaces of 
quasiparticles with different values of rotated-charges [231 ]: 

tt + i = diag(l, 1,0, 0,0,0), (2.4) 

fi+i = diag(0, 0,1, 0,0,0), (2.5) 

Q_i = diag(0, 0,0, 1,1,0), (2.6) 

tt = diag(0, 0,0,0, 0,1). (2.7) 

This is a complete set of projectors, satisfying the following relations: 

QqQq, = 6qq,Qq, Q, Q' = ±1/2, +1, 0, (2.8) 

n + i + n+i + n_i + n = l- (2.9) 

By making use of these projectors, we can decompose the multicomponent quark spinor field 
into separate pieces, describing groups of quasiparticles with different rotated charges: 

ip = +"0(0), (2.10) 

where, by definition, 

= fi+i^j ^(+1) = ^+1^, ^(-i) = fi-i^j V'(o) = Skip- (2.11) 

In the new notation, the quadratic part of the quark Lagrangian density can be rewritten 
as follows: 

= E ^(il^-m + ^ + eQ^A^Qy (2.12) 

Q=±l/2,+l,0 

As follows from Eq. (12.31) . the chemical potentials (jlq for quasiparticles with different Q- 
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charges, when projected onto the relevant color-flavor subspaces, are given by 

2 1 . . 

(2.13) 

(2.14) 

(2.15) 

(2.16) 

In this study, in order to simplify the analysis of the gap equation we will eventually neglect 
the effects due to nonzero /i e and yUg- This is certainly justified in the study of QCD at 
asymptotically large densities. On the other hand, if the analysis is to be extrapolated to 
moderately large densities, relevant for compact stars, nonvanishing fi e and fis m &y become 









= M~ 




3*" 




= Hdr 


= Hdg 


= /i + 
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^8, 


^(+i) 




= fl- 
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2 

3^8, 




M(0) 


= fJ'db 


= M + 


1 


2 

3^8- 





important 
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201 ]. One should keep in mind, however, that the study of such a moderate 



density quark matter from first principles will be still quantitatively unreliable within the 
framework of the Schwinger-Dyson equation because of the strong coupling regime. As for 
the main purpose of this study, it aims only at a better understanding of the qualitative role 
of long-range forces. 



III. QUASIPARTICLE PROPAGATORS 

In the 2SC color superconducting phase, only the quasiparticles with the charges Q = ±~ 
participate in Cooper pairing, while the remaining two quasiparticles (with charges Q = 0,1) 
play the role of passive spectators. Therefore, in the rest of the analysis, we will concentrate 
exclusively on the two pairs of quasiparticles participating in Cooper pairing and ignore the 
others. 

As usual in studies of color superconducting phases, it is convenient to introduce the 
Nambu-Gorkov spinors, 



(-Q) 



for quasiparticles with the charges Q = ±~. Here ib^- = Cibf- and ib^ = ibJ A .C are 

* \Q) (Q) (Q) (Q) 

the charge-conjugate spinors, and C = ry 2 7° is the charge-conjugation matrix satisfying 
the relations: C~ lr y^C = — (7^) T and C = —C T . In terms of the Nambu-Gorkov spinors, 
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Lagrangian density (I2.12p takes the form 



(Q)' 



(3.2) 



Q=±l/2 



-1,0 



where the inverse free propagator S - 1 for each sector with a fixed value of Q-charge has a 
block-diagonal form, 



and the explicit form of the diagonal elements reads 



-i 



(O),o. 



-i 



m. 



(3.3) 



(3.4) 



For quasiparticles participating in Cooper pairing, the full propagators also have nonzero 
off-diagonal Nambu-Gorkov components, determined by the color superconducting gap func- 
tion, i.e., 



rr-l 

(Q) 



[G hJ 1 A (Q) 



(3.5) 



The color-flavor structures of A. A , and A^-. are given by 

(Q) (Q) to J 



-«7 5 A 
27 s A 

i-f 5 A* 
-«7 5 A* 

Note that the explicit forms of the two relevant Nambu-Gorkov spinors (13.11) read 



A 



A 



(+§) 



(+§) 



- A K) 



- A H) 



(3.6) 
(3.7) 



(+1) 



4>ug 



(3.8) 



It appears that one can partially diagonalize the inverse full propagators by simply 



reordering the components of the spinors as follows: 



*<+i> 



i'% 

4>ug 



'rug 
^dg 



(3.9) 
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From physics viewpoint, the possibility of such a partial diagonalization reflects the fact that 
there are two different types of Cooper pairs: one made of red up and green down quarks 
and the other made of green up and red down quarks. 

In the new basis, the inverse full propagator S 1 "- 1 has the following block-diagonal form: 



x 1-1 ro^ -1-1 
(Q) J 



where 

ft 

and 



X i-l 

(+|) J 



( 



\ 



-i 7 5 A* 



i 7 5 A* 



i 7 5 A* 



-i 7 5 A 



7 m (m^ + ~eA M ) + /i u;?7 - m 



i 7 5 A 



i 7 5 A 



^{idf, - \eA^) - jj ug <y° - 



rn 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



7 ^(i(9 M - \eA^) + /i d97 - m 



— i 7 5 A 



-z 7 5 A* 



(3.14) 



Using the representation for the inverse quasiparticle propagator in Eq. (I3.10p . we find the 
propagator itself, 

s (q) = d ^( s 5)> s c$)) • ( 3 - 15 ) 



The calculation of the corresponding diagonal blocks S,^ is tedious, but straightforward. 
The details of derivation are presented in Appendix lAl 



IV. GAP EQUATION 

In the coordinate space, the gap equation (i.e., the off-diagonal component of the 
Schwinger-Dyson equation for the full propagator) reads 



(Q) 



[u, u 
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') = f*V (~T A ) 



aX 

(Q) 



J 21 



(u,u')YT B D* B (u,u'), 



(4.1) 



where D^ u (u,u') is the gluon propagator, and u = (t,z,r±) is a four-vector of space-time 
position. We will assume that the gluon propagator is diagonal in adjoint color indices. Note 
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that the off-diagonal component of the propagator S?q. satisfies a similar equation. While 
Eq. (14. ip describes Cooper pairing of red up and green down quarks, the equation for S 1 ^ 
describes Cooper pairing of green up and red down quarks. 

In this study of Cooper pairing in a magnetized color superconducting phase, it is conve- 
nient to start from the coordinate-space representation of the gap equation [see Eq. ( 14. ip ] and 
then switch to the Landau-level representation. This is in contrast to the usual momentum 
space representation, often utilized in the case of vanishing external fields. 

In this connection, a short remark is in order regarding the general structure of a quasi- 
particle propagator. Because of the interaction of charged quasiparticles with the magnetic 
field, their momenta in the two spatial directions perpendicular to the field are not well- 
defined quantum numbers. This is reflected in the structure of the propagator (as well as its 
inverse), which is not a translationally invariant function in coordinate space. Instead, the 
quasiparticle propagator has the form of a product of the universal Schwinger phase (which 



spoils the translational invariance) and a translationally invariant part .32] (for details, see 
Appendix lAl). 

After factoring out the same Schwinger phase on both sides of the gap equation and 
projecting the resulting equation onto subspaces of different Landau levels, one obtains an 
infinite set of coupled equations; see Eq. (1B6|) in Appendix [B] For both charges Q = ±1/2, 
the gap equations are similar. Here we show only the final set of equations for Q = +1/2: 



WO) ^nj, , WO) 



x YD»u(u-u;',k 3 -k' 3 ;q ± ), (4.2) 

where m, n = 0, 1,2, . . . are Landau-level indices, and functions C n and £ n are defined in 
Appendix [Aj see Eq. (1A19|) and (1A23|) . respectively. These functions depend on the param- 



eters of the model (e.g., masses and chemical potentials of quarks) as well as on the color 
superconducting gap parameters A n . Note that the gaps associated with different Landau 
levels are not necessarily equal. This fact is emphasized by the Landau-level subscript n in 
the notation. Here and below, we assume that all gaps A n are real functions. 
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A. Gluon propagator 



In dense quark matter, unlike in vacuum, the gluon exchange interaction is partially 
screened. Therefore, when analyzing the Cooper pairing dynamics between quarks, it is very 
important to take the relevant screening effects due to nonzero density into consideration 
8) . In the problem at hand, in addition, one should account for the external magnetic field, 
which can further modify the screening of the one-gluon interaction through quark loops. 
The latter can be quite important in strong magnetic fields 33] . To simplify the analysis in 
this study, we will assume that the magnetic field is weak (\eB\ <C fi 2 ). At the end, we shall 
see that this happens to be a very good approximation for most stellar applications. 

In the case of a weak external field, the screening of the one-gluon interaction in dense 
medium can be described well by the usual hard-dense loop approximatio n [34 361 . In the 
Coulomb gauge, the Lorentz structure of the gluon propagator is given by [37, 



Q 2 SfjaS^o P T 



flu 



(4.3) 



q 2 Q 2 -F Q 2 -G" 

where functions F and G define the spectra of the longitudinal and transverse gluons, re- 
spectively. Both functions depend on the energy q° and the absolute value of the three- 
momentum \q\. By definition, Q = (q°,q) is a momentum four-vector. The transverse 



Lorentz projector P^ u is defined as follows: 



P T 



p T — 



P T = 5- 

1 13 u is 



(4.4) 



In the most important regime for Cooper pairin g dy namics, q° <C |g| <C mj), the approximate 



expressions for these screening functions read 

F ~ m 



3J-[36| 



D- 



r< 71 2 Q 
G ~ —m D — 



(4.5) 



where m 2 D = (gfi/ir) 2 is the Debye screening mass in two-flavor quark matter. At large den- 
sities, the exchange interaction by electric gluon modes is strongly suppressed due to Debye 
screening and, to leading order, plays no role. Magnetic gluon modes, on the other hand, 
are subject only to a mild dynamical screening^ (Landau damping) at nonzero frequencies 
and play the dominant role in Cooper pairing 
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B. Gap equation: Weak magnetic field limit 



In order to obtain the gap equation in the weak magnetic field limit, we expand the 
translationally invariant part of the full fermion propagator in powers of the magnetic field 
and keep the leading terms up to second order, (eQB) 2 (for details, see Appendixes O and 
ID]) . Omitting the technical details, here we present the final form of the gap equation, 



A(w) = r<°>(w) + r (1) (w) + t( 2 \u) 



(4.6) 



where 



/At,)' f rl 3 W 



(4.7) 



3 J 2tt J (2tt) 3 

is the contribution of the ith order in powers of the magnetic field. The explicit form of 
the kernels K^\co, k) for the three leading-order terms in the gap equation are presented in 



Eqs. (EH, flC23|) and f lC24j) in Appendix 
At zero magnetic field, Eq. ([ 



6]) reduces to the well-known gap equation in the 2SC phase 



without a magnetic field js-ll4j]. After switching to the Euclidean space and performing 
the traces on both sides of the gap equation, we rederive the following zeroth-order (i.e., 
vanishing magnetic field) equation: 



A (o) ( 



duj' E 

W) 



d 3 k' 



(2tt) 3 ((J ^ 2 



X 



Ll +(#-a0 2 + [a<«K)] 

2|fe-fc'| 
+ \k - k'\ 3 +cu 3 



where k = \k\, k' = \k'\, ue 



l(u E -uj' E ) 2 + \k-k>\ 2 + m 2 D 

icj', and uf = (7r/4)m|,|a;^ 



i J 



(4.8) 



10J, OJ i 



uie\- For simplicity, 



here we assumed that m = and that the chemical potentials of all quarks are identical and 
equal ft. 

After performing the integration over k' and keeping only the leading-order contributions 
from the dynamically screened magnetic gluon exchange, we arrive at 



367T 5 



AW, 



ijj. 



In 



A 



VK) 2 + (AW)2 | W / 



we 



(4.9) 



where A = 4(2/i) 3 /(7rm 2 :) ). The approximate solution to this equation reads 8Ml4| 



A' ' ~ Aexp(— ^- + 1). 



(4.10) 



11 



Using this result as a benchmark, let us proceed to the case of a weak but nonzero magnetic 
field. 

It is easy to check (and might have been expected from the symmetry arguments) that 
the first-order term, i.e., T^'(u) in Eq. ( 14. 7p . which is linear in a magnetic field, vanishes 
after the Dirac traces are performed. Thus, the leading correction to the gap equation in a 
weak magnetic field comes from the second-order term, i.e., T^ 2 >{lo) in Eq. (14. 7p . 

To the same leading order in coupling, which includes only the exchange interaction due 
to dynamically screened magnetic gluons, we derive the following explicit form of the gap 
equation (for the details of derivation, see Appendix iDl) : 



to 



E) 



1 , A 

hi^ r + 



9cul 5 (eQB) 2 sm 2 9 Bk . co t 



In 



(4.11) 

4fi 2 [uf + [«) 2 + (A<*>)= 

The detailed analysis of this equation may not be very easy. However, several of its qualita- 
tive properties are obvious right away. First of all, the positive sign of the subleading-order 
correction, proportional to (eQB) 2 , indicates that the gap increases with the magnetic field. 
This is in qualitative agreement with the intuitive expectation that the external magnetic 



field should en 

charges 
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ranee the binding energy of Cooper pairs made of quasiparticles with opposite 
4231] . From the fact that this correction to the gap equation is also proportional 
to sin 2 #Bfc, where 9sk is the angle between the quasiparticle momentum and the magnetic 
field, we conclude that the gap function acquires a directional dependence. Moreover, we see 
that the largest value of the gap will be for quasiparticles with the momenta perpendicular 
to the magnetic field. On the other hand, for quasiparticles with the momenta parallel to 
the field, there is no enhancement of the gap at all. 

In order to understand the qualitative effect of the subleading term quadratic in magnetic 
field, we can perform the following semirigorous analysis of Eq. (14.111) . To this end, let us 
cut the infrared region of integration off at co' E = A^ B ' and substitute A^ B ' = in the 
denominators of both terms on the right-hand side of the equation. We then arrive at 

18tH y #V J J A (b) u' E u' E 

While this approximation cannot be used to get a reliable estimate for the gap, it is very 
helpful to understand the qualitative effect of the magnetic field on the pairing dynamics in 
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color superconducting dense quark matter. It shows that the effective coupling constant in 
the presence of a magnetic field becomes larger, i.e., 

9 2 -+ & = 9 2 {l + sin 2 B )j , (4.13) 

where we substituted Q = ±|. The validity of the weak field approximation requires that 
the subleading correction is small compared to the leading result. This translates into the 
requirement \eB\ 2 < g 2 ji A . As we shall see below, this condition is always satisfied in stellar 
applications. 

Without rigorously solving the gap equation (14. lip , now we can claim that the solution for 
the gap function in the magnetic 2SC phase is approximately given by the same expression 
as in the absence of the field, but with the coupling constant g replaced by g e g, i.e., 

3tt 2 

A^~Aexp( jL — + 1) ~ A«V Bfc , (4.14) 

where the explicit expression for fish follows from Eq. ( 14.131) . 

fi » = W** f9, » (4 - i5) 

This is a nonnegative function, which depends on the angle between the quasiparticle mo- 
mentum k and the magnetic field B. Its maximum value is obtained at 9bu — 90°. 

The final result in Eq. ( I4.14p is interesting for several reasons. Most importantly, it shows 
that the gap is nonisotropic, taking its largest values when the quasiparticle momenta are 
perpendicular to the direction of the magnetic field, and taking its smallest value when the 
quasiparticle momenta are along/ against the field. We also find that, compared to the case 
without the magnetic field, the gap is subject to an increase in all directions of quasiparticle 
momenta, except for the directions exactly along or against the magnetic field. 



C. Gap equation: Strong magnetic field limit 

To get a qualitative insight about the pairing dynamics in the case of a strong magnetic 
field, eB > /i 2 , it seems sufficient to consider the gap equation in the lowest Landau- level 
approximation. The choice of a simple approximation for the gluon exchange interaction 
is much harder to justify. Here we will use the gluon propagator with the screening effects 
at zero magnetic field. Obviously, such an approximation is not very reliable. A naive 
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justification for such an approximation is the observation that gluons couple not only to 
the charged quasiparticles (with Q = ±1/2 and Q — 1), which are strongly affected by 
the magnetic field, but also to neutral quasiparticles (with Q — 0), which are not affected 
by the magnetic field at all. If zero density (/i = 0) and strong magnetic field limit in 
gauge theories is used as a guide for intuition, one may suggest that those gluons, which 
are coupled only to charged quasiparticles, will be subject to an additional Debye screening 



with an effective mass m e ^ oc gy \eB\ 33]. The other gluons will be still providing the same 
dominant interaction with dynamical screening as in absence of the external field. Then, 
the usual hard dense loop approximation may be still qualitatively reasonable. Besides, to 
the best of our knowledge, the explicit result for the polarization tensor (screening) in dense 
QCD matter (/i ^ 0) in a magnetic field (B ^ 0) is not available in the literature. Thus, 
the main purpose of our exercise in this subsection, which is based on the simplest possible 
approximation, will be to roughly estimate the color superconducting gap due to long-range 
interaction in the regime of a strong external magnetic field. 

By making use of Eq. (14. 2L we easily derive the gap equation in the lowest Landau-level 
approximation, 

du' E dk' 3 f d 2 q ± ( q\l 2 \ & {B) {u' E ) 



{B) g 2 f dco' E dk' 3 fd 2 q ± ( q\l 2 \ 



(co' E ) 2 + (A;' 3 - /x) 2 + (A( B )) 



2 



qj j(k> 3 - k 3 ) 2 + q 2 ^ 

(fc*-fc3)2 + g 2 [{k ,3 _ k 3 )2 + q 2 ±] l + ^ ^ ■ > 

Because of the exponential suppression in the integration over the transverse momentum q±, 
the dominant contribution comes from the region of small momenta, q±l < 1. Therefore, 
an approximate result can be obtained by simply making a sharp ultraviolet cutoff at q± = 
\/2/l and dropping altogether the exponential factor exp(— g 2 _/ 2 /2) in the integrand. After 
performing the integration also over the longitudinal momentum k' 3 , we will arrive at the 
following approximate gap equation: 

727r 2 i_ 00 ~~V(c4) 2 + (A( B )) 2 n K-^l 

where A B = 8 ^, 3 = 87rv/ ^<3-^l 1 . As we see, this equation has the same structure as 
Eq. (14. 9p . but with a smaller effective coupling and a different expression for A#. Making 
use of this fact, we can get an approximate solution for the gap in the limit of strong magnetic 



i m W»L ,..J KU 'L »>rT^7T. (4-17) 
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field by properly modifying the result in Eq. (I4.10p . i.e., 

Af m An\BB\ 3 / 2 , 3vr 2 . . 

A (B) = ' ' exp + 1 . 4.18 

9 V 9 

Here we substituted Q — ±|. This result shows that the strong magnetic field strengthens 
the diquark pair formation. This is in qualitative agreement with the findings in models 



with local interaction 
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In contrast to the result in the weak magnetic field limit, there is no directional depen- 
dence in the gap function when the field is strong. This suggests that the corresponding 
pairing dynamics is essentially local. While the result may appear surprising at first sight, 
this finding in fact agrees with the intuitive picture that the motion of charged particles is 



restricted over distances of the order of the magnetic length, I = 1/ y\eQB\, in the plane 
perpendicular to the magnetic field. When Cooper pairs form, the additional spatial restric- 
tion on particles' motion (partial localization) can strongly enhance the binding energy and 
substantially reduce the size of bound states. 



V. CONCLUSION 



non In this paper, we studied the effect of a rotated magnetic field on the Cooper pairing 
dynamics in the two-flavor color superconducting phase of dense quark matter with long- 
range interaction provided by the one-gluon exchange with dynamical screening. Using the 
Landau-level representation, we derived a set of gap equations valid for an arbitrary magnetic 
field. These equations show that, in general, the gaps are functions of the Landau-level index 
n. Therefore, solving the corresponding set of equations may be rather involved and require 
the use of sophisticated numerical methods. Instead, here we used analytical methods to 
investigate the limiting cases of weak and strong magnetic fields. 

In the weak magnetic field limit, the energy separation between the Landau levels is 
vanishingly small and there is no reason to expect a strong dependence of the gaps on the 
corresponding discrete index n. This justifies the use of an approximation in which the 
gaps are the same in all Landau levels near the Fermi surface. Additionally, in this case 
the quasiparticle propagator allows a simple expansion in powers of the magnetic field that 
greatly simplifies the structure of the resulting gap equation; see Eq. (I4.6p . We find that the 
leading-order term, affecting the gap, is quadratic in the magnetic field. The corresponding 
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correction to the vanishing magnetic field result for the gap is determined by the value of 
parameter /3^ ax ' ) = 81n 3 (eB) 2 / (A\^2g 3 p, A ) , where B is the magnetic field and ft is the quark 
chemical potential; see Eq. (14. 14ft . The numerical value of this parameter appears to be quite 
small even for strongest possible magnetic fields in compact stars, B < 10 18 G. Indeed, the 
corresponding numerical estimate reads 



(Here, for the strong coupling constant, we used g = v4tt, which corresponds to a s = 1.) 

The most interesting feature of the pairing dynamics in the presence of a magnetic field 
is a directional dependence of the gap function in momentum space. The magnetic field 
correction to the gap is proportional to sm 2 9 Bk, where Qbu is the angle between the quasi- 
particle momentum k and the magnetic field B. From the physics viewpoint, this means 
that quasiparticles with momenta pointing perpendicular to the direction of the magnetic 
field have the largest gaps, while quasiparticles with momenta along/ against the field have 
the smallest gaps. Clearly, such a directional dependence is a qualitative outcome of a long- 
range interaction in the model used. This contrasts with the studies based on models with 



pointlike interactions in Refs. |2ll-l23|. where the gaps are always isotropic. 

Our analysis in the case of a strong magnetic field is admittedly less rigorous. We use 
the lowest Landau-level approximation and utilize the simplest approximation for the gluon 
exchange interaction without modifying the screening effects due to a nonzero magnetic 
field. The resulting estimate for the gap is given in Eq. (I4.18p . Our result shows that strong 
magnetic fields enhance the diquark Cooper pairing and lead to larger color superconducting 



gaps. This is in qualitative agreement with the findings in Refs. 2jJ-|23|, where the models 
with short-range interactions were used. We also find that, because of the partial localization 
of quasiparticles in a strong magnetic field, the corresponding dynamics is essentially local 
and there is no directional dependence of the gap. 

To go beyond the two limiting cases, analyzed in this paper, one will need to properly 
truncate an infinite set of gap equations and use numerical methods to solve it. In such 
an approach, it may be also possible to include the effects of different quark masses and 
chemical potentials. The corresponding study, when extrapolated to the regime of realistic 
densities, may further extend our understanding of dense quark matter by clarifying (i) pos- 
sible directional dependences of the gap function, (ii) the evolution of such a dependence 
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between the two limiting cases studied here, and (iii) the effect of /3-equilibrium and neu- 
trality of quark matter on the gap function in magnetic fields. All of these topics are left 
for future investigations. 



Acknowledgments 

The authors would like to thank E. Gorbar for interesting discussions and useful com- 
ments. This work is supported in part by the U.S. National Science Foundation under Grant 
No. PHY-0969844. 



Appendix A: Quark propagator 

In this appendix we calculate the explicit forms of the full propagators for quasiparticles 
with Q = +| charge. (The result can be also easily generalized to quasiparticles with 
Q = — | charge.) We present the details of the analysis for 11- and 21-components of the 
propagator Sf ... 

The starting point of the derivation is the definition of the inverse propagator in 
Eq. (13. lip . Introducing a shorthand notation for the diagonal and off-diagonal elements 
of that propagator, we write 



qX 



A" 



[G 



-l-i 



(Al) 



where 



G ± = [(Gq) - A T G^A 



±1-1 



-G^A ± G ± . 



(A2) 
(A3) 



The explicit forms of the 11- and 21-components of the propagator read 



^ -fi dgl ° + m)-A 2 



-i -i 



qX 

D (+|)21 



-i 7 5 A* 



+ /i„r7 - m j 17 % - Vdgl +m 



A 



-i 



(A4) 
(A5) 
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where, by definition, Trji = id^ + eQA^ and the gauge field is N 1 = (0, 0, xB, 0) with the 
strength of the external (rotated) magnetic field denoted by B. 

The inverse of the operator in the square brackets of Eqs. (IA4I) and (IA5I) . which is the 
same for all components of the propagator, can be calculated by employing the usual trick 
of "quadrating" the operator. In this case, however, we end up "biquadrating" it because 
the corresponding operator is already quadratic in energy For this purpose, let us introduce 
the following shorthand notation: 



X 



± 



(id t - SfiY - TTi - ieQBYr - ) —m — pr — A 
± 2 7 °/i ( 7± • tt ± + 7 3 vr 3 - m) , (A6) 



where Sfi = p = n ± = (n\n 2 ) and 7± = ( 7 \ 7 2 ). Note that X~ is the 

same operator that appears in the square brackets of Eqs. flA4j) and (1A5I) . For simplicity of 
notation, we dropped index Q here. 

Let us first concentrate on the 11-component of the propagator. It can be rewritten as 
follows: 

= (A- 1± -tz ± b)c-\ (A7) 
The three new operator functions introduced here are defined by 

A = [(tdth - vrV - /i dg7 + m\ [{idt - <^) 2 - Ji 2 - 2 / x( 7 3 vr 3 + m) 7 ° - (vr 3 ) 2 - m 2 - A 2 ] 

- [(id t )j° - vr 3 7 3 + /i ur7 ° + m] (tt 2 + zeQ^yV) , (A8) 
£ = (i$ - /i d9 ) 2 - Tri - ieQB^j 2 - (tt 3 ) 2 - m 2 - A 2 , (A9) 
C = - c^) 2 - tt^ - ieQBl 1 ! 2 ~ (tt 3 ) 2 - m 2 + ^ 2 - A 2 j * 

- 4/i 2 [(iSfe - Sfxf - A 2 ] . (A10) 

In the coordinate space, the corresponding propagator is formally given by 

S£ + x )n (u,u') = (u\ (A - 7j. • tt ±j b) C- l \u'), (All) 

where w = (t, z, r±) and r_i_ = (x, y). It is easy to perform a Fourier transform in time and 
^-coordinate, 

Sf + ± )n (u,k 3 ;r ± ,r' ± ) = j ' dtdze iujt ~ ik3z S*_i )u («,«'). (A12) 
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In essence, this transform results in a simple replacement of idt — > uj and tt 3 — > k 3 in all of 
the earlier expressions. 

To proceed further, we should find a basis of suitable eigenstates, in which the propagator 
has the simplest possible form. To this end, we note that the functions A, B and C depend on 
the operator tz\ +ieQB r y lr y 2 . Its eigenvalues are well known: 2n\eQB\, where n = 0, 1, 2, . . . 
is the Landau-level index. Note that the integer quantum number n has both the orbital 
and spin contributions, i.e., n = k + (1 + s)/2, where k — 0, 1, 2, . . . labels a specific orbital 
state, while s = ±1 corresponds to a given (up or down) spin state. The explicit form of 



the corresponding eigenstates (r±\kp y s) is also well known (e.g., see Ref. 39|, where similar 
method and notations are used). 



Following closely the approach of Ref. 



we use the complete set of eigenstates to 



simplify the expression for the propagator flA12j) . The final result will have the form 



Sf +k)21 (u, k 3 ; P jl, r' ± ) = e^ r ^Sf +1 _ )21 (u, k 3 ; r ± - r' ± ), (A13) 

where $(rj_, r' ± ) is the Schwinger phase. In the Landau gauge used, the explicit form of the 
phase is 

*(ri, ri) = - (X + ^~ ?//) sgn(gg^), (A14) 



where / = 1/ y \eQB\ is the magnetic length. (Note that this phase is responsible for 
breaking the translational invariance of the propagator.) The translationally invariant part 
of the propagator is given by 

Sf H)11 (co,k 3 ;r ± ) = ^ { ^ + 2^_i(0^+] - i^^^^JUi^) j fA15) 

where ^ = r 2 L /(2/ 2 ), are the generalized Laguerre polynomials (by definition, L n = L° 

and L" 1 = 0), and 

V± = i (l ± « 7 1 7 2 sgn(gg5)) (A16) 

are the spin projection operators. 

Functions A n , B n and C n in Eq. (|A15[) replace the corresponding operators A, B and C, 
when projected onto the nth Landau-level state. Their explicit forms are obtained from A, 



19 



B and C by replacing irj_ + ieQB r y lr y 2 — > 2n\eQB\, i.e., 

An = [ojj° - fc 3 7 3 - W ° + m] [(to - 5 fx) 2 -ft 2 - 2p{ 7 3 k 3 + m) 7 ° - (A; 3 ) 2 - m 2 - A 2 ] 

- 2n\eQB\ [^7° - fcV + /i ur 7° + m] , (A17) 
B n = (w - /i dg ) 2 - 2n\eQB\ - (A; 3 ) 2 - m 2 - A 2 , (A18) 
C n = \(u - 5ll) 2 - 2n\eQB\ - (k 3 ) 2 - m 2 + ft 2 - A 2 1 ' - 4/i 2 [(to - Sli) 2 - A 2 ] . (A19) 

Here we consider a general case when the dynamically generated gap function A n depends 
not only on the energy u and k 3 , but also on the Landau- level index n. (In operator form, 
it means that A depends on tz\ + ieQBj 1 '-/ 2 .) Therefore, we replaced the operator A with 
the corresponding value A„ that it takes in the nth Landau-level state. 

At this point, it may be appropriate to note that the zeros of C n determine the spectrum 
of quasiparticles in color superconducting quark matter in a magnetic field, i.e., 



E n ,±,± =Sfl± 



2n\eQB\ + (A; 3 ) 2 + m 2 ± p, 



A 2 . (A20) 



Note that all four different sign combinations are possible. The choice of the sign in front of 
the chemical potential p, corresponds to the choice of either particle states (allowing small 
energies of order A n ) or antiparticle states (generally having large energies of order p). The 
sign in front of the overall square root corresponds to particle/hole type quasiparticles (i.e., 
positive /negative energy states). One should note, however, that an additional complication 



in 



;his classification appears in the case of gapless superconducting phases when 5fi > A r 
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19] . 



Following the same approach, we can derive explicit expressions for all components of 
the propagator Sf_ [ iv For example, the final expression for the off-diagonal 21-component, 
which is used in the gap equation in the main text, reads 

Sf +h)SL (u, k 3 ; r ± , r' ± ) = e^ r ^ Sf +1 _ )21 (oo, k 3 ; r ± - r' ± ), (A21) 

with the translationally invariant part given by 

—£/2 00 fa 9— ^ 

(A22) 



Here we introduced yet another function, 

S n = {u- 5fi) 2 - 2n\eQB\ - (k 3 ) 2 - m 2 - p 2 - A 2 - 2p{k 3 -f 3 + 771)7°. (A23) 
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Before concluding this appendix, let us add that similar representations can be also derived 
for the components of the inverse propagator. As an example, let us present the correspond- 



ing result for 



qX 



21 



(u,u f ), which is used in the gap equation. It has the same general 



structure as the above expressions for the components of Sf Lv i.e., 



r-A 



2>J 21 



u,k 3 ;r ± ,r' ± 



i<5>{r ± ,r' ± ) 



S x 



21 



u, k 3 ; r ± - r',). 



(A24) 



It is important that the inverse propagator has exactly the same phase as the propagator 
itself; see Eqs. (1A13|) and f ]A14j) . The explicit form of its translationally invariant part reads 

1 p-C/2 00 



s x 



(+§) 



21 



-27" 



2tt/ 2 



K + 



(A25) 



n=0 



Appendix B: Gap equation 

The gap equation (i.e., the off-diagonal component of the Schwinger- Dyson equation for 
the full propagator) in the coordinate space reads 



9 X ~ 
(Q) 



21 



-w 2 y (t a y 



(Q) 



21 



(u, u')YT B D^(u - u) 



(Bl) 



where D^(u, u') is the gluon propagator, which is assumed to be diagonal in adjoint color 
indices (A, B = 1,2, ...,8), i.e., D^{u — u') = 8 AB D^iu — u'). By making use of the identity 



ErpA rpA 
1 a'a 1 b'b 



A=l 



1 _ 1 
-^a'b$ab' — Q^aa'^bb' 



(B2) 



we derive the following form of the gap equation: 

- 1 , „ .2 



qX 

(Q) 



21 



(Q) 



J 21 



(u,u)^ u D^ v (u - u). 



(B3) 



Taking into account that all components of the quasiparticle propagator as well as its inverse 
have the same nonzero Schwinger phase, we can derive the equation for the translationally 
invariant parts simply by dropping the common phase factor on both side of the gap equation, 
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i , 3 . 2g 2 fdu'dk' 3 - x . . . l 



3 J (2tt) 2 



(2tt)' 



D^(u -u , k - k , q±), 



(B4) 



where we additionally performed a Fourier transform in time and ^-coordinate on both sides 
of the equation, and used a momentum representation for the gluon propagator. 
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By making use of the explicit form of the relevant translationally invariant parts of the 
propagators in Eqs. (IA22j) and flA25j) . we rewrite the last form of the gap equation as follows: 



-C/2 



2tt/ 2 



.2g*e 



2 P -m 



E 



n=0 



dco'dk' 3 „A r 



3 2tt/ 2 ^ 7 (2vr) 2 " 7 C r , 

n=0 ^ V 7 

,-.7± • r ± o 



x 



X 



d 2 q± 



-D^u -u',k 6 - k' 3 ,q ± ). 



(B5) 



The last equation can now be easily projected onto different orbital eigenstates. This is for- 
mally done by multiplying both sides of the equation by e~^' 2 L m (£) (where m = 0, 1, 2, . . .) 
and integrating over the perpendicular spatial coordinates r±. After performing such pro- 
jections, we arrive at the following (infinite) set of gap equations in the Landau-level repre- 
sentation: 



; V ~ f Cbj'dk' 3 f d 2 q ± ^An£n 



3 (2vr) 2 J (2tt)* 



r 



r(o) 



YD^(LO-u',k 3 -k n \q ± 



./3 



(B6) 



where, by definition, 



iy+m e -x L m-n ^ £r 



(B7) 



In the derivation, we used the following table integrals (see formulas 7.414 3 and 7.422 2 in 



Ref. 



4Qj): 



d.nr x x a L a m {x)L a n (x) = r(n + ? + 1) C 



nl 



(B8) 



and 



dxa^ +1 e-^L^(aa?)L*(aa?)J (xy) 



\m+n 



2a 



(7+1 



m + ay. 

■ j — « 

ml 



jn—ml V \ j m—n ( V 



4a 
(B9) 



Appendix C: Propagator in weak magnetic field limit 

In this Appendix, we consider the quasiparticle propagator and the gap equation in the 
limit of weak magnetic field. 
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We begin by performing a Fourier transform of the translation invariant part of the 
propagator, 

Sf +h)21 (co,k 3 ,k ± ) = f d 2 r ± e- ik ^S x + i )21 (u,k 3 ;r ± ) 



00 AC 

-2z 7 V^ 2 £(-l) B =^[L». (2fciZ 2 ) P_ - L n ^ (2kll 2 ) V, 

n=0 v 

+4^-^)7°^ (2A^ 2 )}. 



(CI) 

In the weak field limit, the difference between the neighboring levels is vanishingly small in 
energy and the properties of the corresponding states become almost indistinguishable. In 
application to the gap function A n , this means that it will become almost independent of 
the Landau-level index in a wide range of n near the (would be) Fermi surface. (Strictly 
speaking, the true Fermi surface is not well defined in a superconductor, but if the gap is 
small, A <C p,, one could map the corresponding phase space onto the phase space in the 
free quark matter.) 

In order to derive a weak field expression for the propagator, one needs to first perform 
the sum over the Landau-level index n. A straightforward way of achieving this is to employ 
the usual proper-time representation, i.e., 

cw^ 'ft^' (C2) 

= fcfacosMe-^H) (03) 

for the two types of structures appearing in the Euclidian propagator, and then use the well 
known summation formula for Laguerre polynomials, 

OO , V 

K(x> n = (1 - z)- [a+1) exp • (C4) 

Before using these identities, it is convenient to rewrite propagator (1C1I) in the following 
form: 

Sf +1)21 (iu E , k 3 , k ± ) = * 7 5 A [I x + 2/i(A; 3 7 3 + m + m°h°h + 2//( 7 ± • fc±h%] , (C5) 



23 



where, by definition, the sums U (i = 1, 2, 3) are 

Jl " 26 L " I Wj [a + 2(n + l)|6|P + c ^ + ) ( ' C6) 



n=0 

fc2 °° X 2A; 2 \ / 1 1 



n=0 

t 2 oo 



^4 e -wg(- ir4 (^) [o + 2(n + 1)|6|12 + c2 . (08, 



Here we used the following notation: 



a = (co E + i5^f + (A; 3 ) 2 + m 2 + A 2 - /i 2 , (C9) 
6 = eQB, (CIO) 



c = 2/2^/ (uj e + iSfi) 2 + A 2 . (Cll) 

It is appropriate to mention that the use of the proper-time representations, as given by 
Eqs. (1C2j) and (1C3I) . may not be completely justified in the presence of a nonzero density. 
Indeed, when the chemical potential is sufficiently large, the above expression for the param- 
eter a may become negative. When this occurs, the proper-time integrals become divergent 
and the validity of the derivation seemingly fails. The way around this problem is to assume 
that the chemical potential is sufficiently small at all intermediate stages of derivation. In the 
end, after magnetic field expansion is done and all proper-time integrations are performed, 
one can extend the validity of the propagators to large values of the chemical potential. 

With the above remark kept in mind, we use the proper-time representations to rewrite 
the expressions for the sums Ij as follows: 

00 /op \ r°° 

h = 2e -*i/l6l Y(-iyL n 4rf / ds cos(sc)e-< a+2n ^ (V- + e" 2 ^) , (C12) 

IS v i 6 ' / J ° 

j 2 = 2e -ki/\b\ J2(-l) n L n ^ ^ S m(sc)e-< a+2n ^ (P_ + e" 2 ^) , (C13) 

n=0 \ \ \ / JO C 

h = 4e- fc '/l b l fy- l T L n (jjjf) J™ ^ sm(sc)e-^ +2n ^ +2 ^. (C14) 
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Then, after using the summation formula (1C4|) . we derive 

POO 

h = I dscos(sc)e- sa - {k y b)tanHsb) [1 - z 7 V tanh(s&)] , (C15) 
Jo 

f°° rlQ 

h = / -sin( SC )e- sa -«^ tanh ^ [1 -z 7 Vtanh(s6)] , (C16) 
Jo c 



CO 



— sm(sc)e- sa -^ /b)tanh{sb) \—. (C17) 

c cosh (sb) 

Finally, expanding the integrands in powers of the magnetic field b and integrating over the 
proper time, we obtain 

a + k 2 . 



h ~ J ds cos(sc)e- s(a+k ^ (l - z T 1 7 3 s6 + j k\b 2 + 0(b 3 )^j 



[a + k]) 2 + c 2 

[(q + £r[) 2 + c 2 ] 2 [(a + fc 2 ) 2 + c 2 ] 4 



j 2 ~ ^ s in( SC )e-^ a+fc i) (l - i 7 V s& + y A;i6 2 + 0(6 3 )) 

2(a + fc 2 ) 8(a + fc 2 )[(a + fc 2 ) 2 - 

[(a + k 2 ± ) 2 + c 2 ] 2 [(a + k 2 ± ) 2 + c 2 ] 4 

j 3 ~ ^ sin( S c)e- s ( ffl+fc i) - s 2 6 2 + ^ A; 2 ^ 2 + 0(6 3 )^) 



a 



+ /c 2 ) 2 + c 2 



x 2(a + fci) 8(q + fc 2 )[(q + fc 2 ) 2 -c 2 ]fc 2 

1 



(a + fc 2 ) 2 + c 



2 



2[3(a + fc 2 ) 2 -c 2 ] 8(q + fc 2 )[(q + fc 2 ) 2 -c 2 ]fc 2 
~ WTW+W + [(a + A; 2 ) 2 + c 2 ] 4 6 (C20) 

Now, combining the same order terms in powers of the magnetic field, we rewrite propagator 
(1C5P as follows: 

S£i )2 i(«^, A: 3 , k ± ) = z 7 5 A [^ (0) + + AT^] , (C21) 



where 



tn\ at + di, + 4u(7 ■ fc + m)7° , 
2^(0) = fc I k Z ^ LL f (C22) 



2a k a k 

" 77 *' (C3) 
(2) = aiiaj - 4^ 2 ) 3 + (q^ - 4^x 2 ) 3 + Aji 2 a+al [iQjj 2 - 3(q+ + aj] 

- 4/i(7.i ■ fc ± ) 7 ° 4 f 7 !\ 3 a H 2 + 160a fc ( 7 ■ * + m) 7 ° 2 f ~ !* 4 °* k\b 2 . (C24) 
« a k ) « a fc ) 

Note the shorthand notation used, 

7 ■ fc = 7± • fc ± + fc 3 7 3 , (C25) 

at = [oj e + iS/j,) 2 + (E k ± p,) 2 + A 2 , (C26) 

a k = a + k]_ = (ou E + iS/j) 2 + k 2 + m 2 + A 2 - p 2 , (C27) 



25 



as well as E k = \fk 2 + m 2 and k 2 = k\ + (/c 3 ) 2 . 



Appendix D: Gap equation in weak magnetic field limit 

To leading order (i.e., the limit of vanishing magnetic field), the gap equation reads 

A<°>M = jJ^J ^ A(0) K)tr [yjr<°V, *07"] ^ - a/, fc - *0. (Dl) 

Here we assumed that the gap is an explicit function of the energy, but not of the momentum. 
The result for the trace in the integrand is given by 

tr [ 7 ^°V, fc'hl = 2<T a fc , + + 5' + • • ■ , (D2) 

^k'^k' 

where the ellipsis stands for antisymmetric terms, which do not affect the form of the gap 
equation. Indeed, when contracted with the gluon propagator, which is symmetric in Lorentz 
indices, all antisymmetric terms will vanish. 

At asymptotic densities, we can also neglect all corrections due to nonzero m and 5fi. By 
taking into account that the main contribution to the momentum integral on the right-hand 



side of the gap equation comes from the vicinity of the Fermi surface (k' ~ kp = \Jp? — m 2 ), 
we can make the following approximation for the trace: 

tr \YK i0 \u', fc'hl ~ (D3) 

a k> 

Note that, in the vicinity of the Fermi surface, one has 

H> = K) 2 + & + A 2 « (D4) 
a+ = Aji 2 + iflCk' + H, ~ 4/i(/i + Ck'), (D5) 

where = — fi ~ k' — kp. 

The resulting equation coincides with the known form of the gap equation in the case of 
zero magnetic field studied in Refs. js i^]. In our notation, the corresponding solution for 
the gap function reads 

|A (0) | ~ Aexp(-^U + 1), (D6) 
V2g 

where A = 4(2fi) 3 /(irm 2 D ). 
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In order to find the correction to the gap function due to nonzero magnetic field, let us 
include the approximate kernel up to second order in the magnetic field. After taking traces 
on the both sides of the equation, we obtain 

+ jj d -^J ^A( B )(^)tr[ 7 ^ 2 )(^fe0 7 1^(a;-^fc-fe').(D7) 

In addition to the result in Eq. (1D2I) . this equation also contains the trace of the second 
order correction to the kernel. The corresponding approximate expression in the vicinity of 
the Fermi surface reads 

tr [ 7 ^ 2 V, fc')7l * 9^^r^(eQBy + • • • , (D8) 

where Ny — 4/i^/(2^/ — ay) — 24^, + 16a^,^|, — (a^) 2 and the ellipsis denotes antisymmetric 
terms. Let us point out that the only directional dependence of this trace comes through 
the overall factor (k^) 2 = (k') 2 (l — cos 2 dsk'), where Qbw denotes the angle between the 
direction of the magnetic B and the momentum k'. (Strictly speaking, in a self-consistent 
analysis, the gap function on the right-hand side will also have a directional dependence and 
will affect the angular integration. The corresponding effects are expected to be very small 
and will be neglected in the simplified analysis here.) The integrand on the right-hand side 
of Eq. ( TD7j) has an additional directional dependence in the gluon propagator [see Eq. (14. 3p ], 
which is a function of the polar angle 9 = 9kk' (i-e., the polar angular coordinate of vector 
k' measured from the direction of the external vector k). With this convention for angular 
coordinates, it is convenient to use the following relation: 

cos 0Bk' = sin 9 sin 9sk cos(0 — 4>Bk) + cos 9 cos 9bu, (D9) 

in order to rewrite the expression for (k' ± ) 2 in terms of the angular integration variables 9 
(polar angle) and <fr (azimuthal angle). Now we can easily perform the angular integration on 
the right-hand side of the gap equation. The results for the two types of angular integrations, 
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namely with the electric and magnetic part of the gluon propagator, read 
(1 — cos 2 0Bk') sin^ 



7T 



M 2 — 2k'k cos 9 



8(k'fk 3 



1 



2k' kM 2 [1 + 3 cos(2fl Sfc )] 
M 2 + 2k'k 



+ - {A(k'fk 2 [3 + cos(2# Bfe )] — M 4 [1 + 3 cos(20 Bfe )]) In 



A 



mag 



1 — cos 2 0Bk') sin 6 dOdcf) 



2tt 
Wk 



M 2 - 2k'k 

2[(k') 2 + k 2 -2k'kcos6} 1/2 
[(k') 2 + k 2 -2k'kcos6f /2 + u? 

(k' + kf + co? 



(D10) 



M c«,s 2 ^, - ( (fc T,^ 2 V(l-3cos 2 ^.) 



2k'k 



In 



to, 



x 6 dx 



x s + 



T " 2 [(A:? + k 2 ] 



k'-k\ 3 +LOf 

max x 4 <ix 



+ 1 



(Dll) 



where M 2 = (u' E — ooe) 2 + {k') 2 + k 2 + m 2 -,. In order to simplify the calculation of A mag , 
it is convenient to change the integration variable 9 to the new dimensionless variable x = 
(1/uji)sJ (k') 2 + k 2 — 2k'k cos 6. Note that sin6d6 = u?xdx/(k'k) and the new range of 
integration is from £ min = \k' — k\/ui to :r max = + 

In the vicinity of the Fermi surface, the approximate results for these integrals read 

7rsin 2 # Bfcl (2/i) 2 



icl 



In 



/// 



3 



^ ^sin 2 ^ 
mag " 3^ 2 \k>-k\* + uj? 



- u E ) 2 + (A/ - A;) 2 + m 2 , 
(2/i) 



+ 



+ 



(D12) 
(D13) 



where the ellipses denote the subleading terms. 

By making use of the above intermediate results, we arrive at the following form of the 
gap equation, 



1 £J 9 7^ (2tt) 7 (2tt 



In 



(2/i) 5 



fc' - /c| 3 + 



toy 



x 



\ [-24g + 16a^| - (a fc ,) 2 ] (e~Q£) 2 . 
l + ^fe? Sm ^ 



(D14) 



Recall that cof = (ir j 'A)m 2 D \u' E — oje\- Integrating over the momentum, we arrive at 



36tt 2 



E) 



In 



A 



VK) 2 + (A(^)) 2 W E - u E 



+ 



9ul 5 (eQBfsm 2 6 Bk 



Afl 2 ( W f + H)H(AW)f)' 



In 



We ~ u e\ 



(D15) 



28 



To get a rough estimate, let us take an infrared cutoff in the energy integration at 

and drop the dependence on A^ B ) in the denominator of the integrand. Then, we have 

18tt 2 7 a(s) £ K| ^ ttJj, TTij-) J \u' E \ 1 ; 

This means that the magnetic field correction is equivalent to an effective increase of the 
coupling constant, i.e., 

^# c fr = # \ l + -fjjz sm °Bk I , (D17) 

where we used the definition of the Debye mass m 2 D = (gp,/n) 2 . 
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